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We study SU{2) gluodynamics at finite temperature near the deconfining phase transition. We 
create the lattice ensembles using the tadpole improved Liischer-Weisz action. The overlap Dirac 
operator is used to determine the following three aspects of vacuum structure: (0 The topologi- 
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fermion spectral density is determined and found to depend on the Polyakov loop above the phase 
transition and (///) the corresponding locaUzation properties of low lying eigenmodes are investi- 
gated. 
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1. Introduction 

More than ten years ago, using a model generalizing random matrix theory, M. A. Stephanov [jl]] 
predicted, that in SU (3) gluodynamics above Tc the different Polyakov loop sectors behave differ- 
ently. In the complex-valued Polyakov loop sectors the chiral condensate should turn to zero at T 
substantially above T^. For SU (2) lattice gluodynamics, where the Polyakov loop is real, it was 
predicted that the chiral condensate stays non-zero, < 0, for all temperatures T > in the 

sector with a negative averaged Polyakov loop L < 0. 

As for SU (3), Gattringer et al. [Q] came to a different conclusion. They defined a new observ- 
able, the gap in the Dirac spectrum, and used it as an order parameter for the restoration of chiral 
symmetry. It was found that the spectral gap opens up at one single temperature T = Tc in all three 
Z3 sectors. Here we examine whether Stephanov's prediction for the Dirac spectrum remains valid 
in the case of SU (2) gluodynamics in the deconfined phase. 

2. Improved action 

Ensembles of O(IOO) statistically independent quenched SU{2) configurations are generated 
with the tadpole improved Liischer-Weisz action on 20^ x 6 lattices. This action is known to sup- 
press dislocations. The form of the action is : 

S = pimp ^ Spi — 2 51 ' i^-^) 

pi ^^"0 rt 

where Spi and 5^ denote the plaquette and 1^:2 rectangular loop terms in the action. Spirt = 
\Tr{\ — Upirt)- The factor uq = (Wivi)'^^ is the input tadpole factor. It is determined from 
Wijti = {{l/2)TrUpi) computed at zero temperature [j3|]. The deconfining phase transition occurs at 
Pi^p = jSc = 3.248(2) for A^,: = 6, which corresponds to T,■/^/a = 0.71 (2) [§]. 

3. Massless overlap Dirac operator 

The massless overlap Dirac operator has the form ^ 

Dov = ^(^l+Dw/ \/dIDw^ , (3.1) 

where Dw = M — p/a is the Wilson Dirac operator with a negative mass term, M is the Wilson 
hopping term, a is the lattice spacing. The optimal value of the p parameter is found to be 1 .4 also 
for the lattice ensembles under investigation. Anti-periodic (periodic) boundary conditions in time 
(space) directions are imposed to the fermionic field. 
In order to compute the sign function 

Dw I ^D\„Dw = 75 sgn(//iv) , (3.2) 

where H\f = 75 Dw is the hermitian Wilson Dirac operator, we use the minmax polynomial approx- 
imation. The overlap Dirac operator constructed this way preserves the chiral symmetry even on the 
lattice and allows to study the properties of the Dirac modes from first principles. It will be called 
D in the following and replaces the continuum Dirac operator D = 7^ where = — igA^ 
is the covariant partial derivative with the gauge field background A^. 
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Figure 1: Probability distributions of the topological charge Q for two temperatures below and above 7^ . 



4. Topological susceptibility Xtop{T) 



We solved the Dirac equation numerically for its eigensystem 



(4.1) 



and considered 50 lowest eigenvectors. As a first application we search for the exact zero modes. 
Their number is related to the total topological charge Qtop of the lattice configuration through the 
Atiyah-Singer index theorem : 

Qtop = Qindex=N--N+, (4.2) 

where and N+ are the numbers of fermionic modes with negative and positive chirality 75 V'^. 
respectively. For the lattice ensembles the expectation value {Qtop) should vanish, but {Q^op) mea- 
sures the strength of global topological fluctuations. The topological susceptibility is 



Xtop 



topi 



V 



(4.3) 



where V is the four-dimensional lattice volume in physical units. In Fig. [T] (left) we see a histogram 
of the topological charge in the confinement phase, close to the transition. Fig. |I| (right) shows the 
corresponding histogram for a temperature higher up in the deconfinement phase. Both histograms 
can be approximately fitted by Gaussian distributions. 

Let us now discuss the topological susceptibility as function of temperature. In Fig. ^ (left) we 
show, that the topological susceptibility Xtop in the negative Polyakov loop sector (L < 0) agrees 
at all T within two standard deviations with x,op in the positive Polyakov loop sector (L > 0). 
In Fig. ^ (right) we compare our final data for Xtop{T), which are for T < averaged over all 
configurations and for T > Tc only over the subsample with L > 0, with the results of AUes et 
al. [^. These authors presented the values of 10~^ x XtopI I^l various values of j8 for Wilson's 
action representing different temperatures. We took A/, = 14.15(42) MeV and extracted their 
susceptibility Xtop{T) from these data. The topological susceptibility is slowly decreasing with 
increasing temperature for both sets of data. Notice that the overlap definition of Q results in a 
systematically higher susceptibility than the improved field theoretic definition employed by the 
Pisa group. 
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Figure 2: The topological susceptibility Xtop as function of T separately for L > and L < (left), and 
comparison of our final result with that of Ref. (right). 
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Figure 3: The spectral density of eigenmodes of the overlap Dirac operator for two temperatures T < on 
the 20"* X 6 lattice. 



5. Spectral density, chiral symmetry restoration and different Za^ sectors 

The chiral condensate (yy) is related to the density p(A) of the non-zero eigenvalues A at 
A — > via the Banks-Casher relation: 

(mm) = - lim lim ERi^ . (5.1) 

The non-zero modes are globally non-chiral, but the near-zero ones are still locally chiral and 
correlated with lumps of the topological charge density. The number of modes belonging to this 
near-zero band is proportional to the total volume V. In the chirally broken phase the required limit 
(5J) of p(A) is non- vanishing at A = In the chirally symmetric phase one expects p(A) = 



in a finite region around the origin, i.e. that the spectrum develops a gap. For the confinement 
(chirally broken) phase we find indeed that the spectral density in physical units is practically 
constant (almost T independent), as can be seen in Fig. |[ Comparing results for configurations 
with L > and L < we found that at low A the density p (A) is 50 ~ 70 MeV higher for negative 
Polyakov loop sector. We believe that this difference disappears in the thermodynamic limit. 
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Figure 4: The spectral density of eigenmodes of the overlap Dirac operator for four temperatures T > Tc on 
a 20^ X 6 lattice, evaluated separately according to the sign of the averaged Polyakov loop. 
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Figure 5: The spectral gap for SU (2) lattice gauge theory as function of temperature, evaluated separately 
according to the sign of the averaged Polyakov loop. 



For the deconfinement phase, when we take only configurations with an average Polyakov loop 
L > 0, Fig. ^ (left) shows that p(A) gradually decreases with increasing temperature, indicating 
the decrease of the chiral condensate until a gap opens and gets wider. For configurations with 
L < Fig. ^ (right) shows that p(A) at low lambda is nonzero and even grows with increasing 
temperature. 

6. Spectral gap 

The spectral gap was defined by the smallest eigenvalue, which does not belong to a zero- 
mode. In Ref. [Q] Gattringer et al. have shown for SU (3) gluodynamics that the gap, as function of 
temperature, has a similar behavior for the real and both complex sectors corresponding to the phase 
of the averaged Polyakov loop. The phase transition occurs at the same Tc, and with increasing 
lattice volume the gap is decreasing. Analogously, SU (2) gluodynamics has only two sectors in the 
deconfinement phase, distinguished by the sign of the (real-valued) averaged Polyakov loop. We 
show in Fig. |5] a clearly defined and rapidly growing gap for configurations with L > 0, whereas for 
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Figure 6: The average IPR within spectral bins for one temperature T <Tc. 
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Figure 7: The average IPR within spectral bins for four temperatures T > T^, separated according to the 
sign of the averaged Polyakov loop L. 



configurations with L < the gap remains very small up to temperatures several times higher than 
Tc. The small gap is a finite- volume effect and can be made to vanish in the limit of spatial Vt, — > oo. 



7. Localization in different parts of the spectrum 

The scalar density of an eigenmode corresponding to an eigenvalue A is denoted as 

Px{x) = such that LyPa(^) = 1 by virtue of normalization. The inverse participation 

ratio (IPR) Ii is the natural measure of the localization. For any finite volume V it is defined by 

h=Vl^pl{x). (7.1) 

The IPR characterizes the inverse volume fraction of sites forming the support of Px{x). From 
Fig. |6| we conclude that for the temperature near but below Tc the IPR (localization) monotonously 
increases with decreasing eigenvalue. There is no clear mobility edge. ^ The monotony is not 
perfect among the lowest one or two bins. Thus, out of the low lying modes, the higher ones are 
continuously less localized. We found that at these temperatures for configurations with negative 

'For r = the localization of overlap eigenmodes has been investigated in Refs. 
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Polyakov loop, L < 0, the modes are by factor 2 ~ 3 less localized. Again, as for p(A), this 
difference should disappear in the thermodynamic limit. In Fig. we show that with increasing 
temperature the average IPR within the respective eigenvalue bins is increasing. The effect sets 
in for higher and higher eigenvalues corresponding to the mobility edge moving outward in the 
deconfinement phase with the gap for L > 0. In the negative Polyakov loop sector the IPR is 
constant at a low level, except for A < 500 MeV, where the tendency of the IPR to grow exists but 
is very weak. 

8. Summary 

We performed first measurements of the topological susceptibility with the help of the overlap 
Dirac operator in finite temperature SU (2) gluodynamics. We found that the topological suscep- 
tibility in the confinement phase is almost constant and is slowly decreasing in the deconfinement 
phase, in agreement with previous results We did not find systematic effects of the sign of 
the averaged Polyakov loop on the topological susceptibility. The chiral condensate, however, be- 
haves completely different in the L < sector. Chiral symmetry remains broken, the spectral gap 
stays close to zero for all T > Tc in agreement with Stephanov's model predictions. A microscopic 
explanation in terms of the interplay of holonomy and topology needs to be worked of. This 
difference is accompanied by a different localization behavior of the lowest fermionic eigenmodes 
in the two sectors. 
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